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effort  ground  based  measorements  of  OH  airglow  rotational 
temperatures  at  85  km  were  made  by  means  of  an  infrared 
;  Fourier  spectrometer  (IRFWI  of  USU).  These  measurements 

j  employed  a  least  squares  spectral  fitting  technique  involving 

the  3-1  Meinel  band.  These  data  and  their  estimated  PSD’s 
will  be  presented.  The  latter  will  be  interpreted  by  means  of 
the  relations  between  gravity  wave  temperature  fluctuations 
and  vertical  displacement  fluctuations  given  by  Makhlouf  et  al 
(1990). 

The  results  will  be  compared  with  the  predictions  of  the 
pavity  wave  model  of  Dewan  ('990, 1991)  and  shown  to  be 
in  reasonable  agreement  thus  lending  further  support  to  the 
local  wave-cascade  hypothesis. 


i  !wb.  DISTaI3'j  I  iOS 


8 


DTIC 

lELECTE  I 
nECOBISSBI 

^  a  I 


14.  SUBJECT  TERMS 

Gravity  waves,  Airglow,  Middle  atmosphere.  Spectra  wave 


15.  NUMBER  OF  PAGES 

4 


16.  PRICE  CODE 


17  SECURITY  CLASSIFICATION  118.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFir  • -^lON  20.  LIMITATION  OF  ABSTRAC . 
OF  report  of  this  PAGE  OF  ABSTRACT 


Unclassified 


NSN  7540-0 1-280-5500 


SAR 


Sta^carc  ^orn-\  J99  '^ev  C  29 

.-•.ov.  .  ;  -  ,  >%,  A-)'.  '13  ^ 


Unclassified 


Unclassified 


PL-TR-92-2322 


GEOPHYSICAL  RESEARCH  LETTERS,  VOL.  19,  NO.  6,  PAGES  597-600,  MARCH  20,  1992 


MESOSPHERIC  OH  AIRGLOW  TEMPERATURE  FLUCTUATIONS: 
A  SPECTRAL  ANALYSIS 

E.M.  Dewan',  W.  Penclleton%  N.  Grossbard’,  and  P.  Espy^ 


Abstract.  A  field  campaign  under  Project  MAPSTAR  was 
conducted  in  Colorado  during  May-July  1988.  As  part  of  this 
effort  ground  based  measurements  of  OH  airglow  rotational 
temperatures  at  85  km  were  made  by  means  of  an  infrared 
Fourier  spectrometer  (IRFWI  of  USU).  These  measurements 
employed  a  least  squares  spectral  fitting  technique  involving 
the  3-1  Meinel  band.  These  data  and  their  estimated  PSD’s 
will  be  presented.  The  latter  will  be  interpreted  by  means  of 
the  relations  between  gravity  wave  temperature  fluctuations 
and  vertical  displacement  fluctuations  given  by  Makhlouf  et  al 
(1990). 

The  results  will  be  compared  with  the  predictions  of  the 
gravity  wave  mode'  of  Dewan  (1990,  1991)  and  shown  to  oe 
in  reasonable  agreement  thus  lending  further  support  to  the 
local  wave-cascade  hypothesis. 

Introtluction 

A  field  campaign  was  conducted  in  Colorado  during  May 
through  July  1988  under  project  MAPSTAR  of  AFOSR.  Its 
purpose  in  part  was  to  study  the  effects  of  gravity  waves  upon 
the  temperature  of  the  OH  airglow  layer  located  around  85  km 
altitude.  To  do  this  we  employed  the  IRFWI  Fourier 
spectrometer  of  Utalt  State  University  from  which,  by  means 
of  the  3-1  vibration  transition  rotational  lines  in  the  Meinel 
band  (near  1.53  p),  we  obtained  the  OH  rotational 
temperatures  by  employing  a  technique  similar  to  that 
described  in  Hill  et  al.  (1979).  The  instrument  had  a  field  of 
view  of  0.8°  and  a  resolution  of  2.5  wave  numbers.  Due  to  the 
high  throughput  of  the  IRFWI  instrument,  temperatures  of  the 
OH  Meinel  band  in  question  could  be  extracted  at  around  a  30 
second  temporal  resolution  with  a  relative  error  of  between 
+  0.5°  K  and  +  1.5°  K  for  each  of  the  23  nights  of 
observation.  For  pointing  angles  both  zenith  and,  more  often, 
25°  elevation  angles  were  employed. 

The  main  purposes  of  the  present  paper  are  to  exhibit  the 
power  spectral  densities  (PSD's)  of  these  temperature  fluc¬ 
tuations  and  to  compare  them  with  the  theoretical  expectations 
based  on  Dewan  ( 1 990,  1991). 

Temperanue  Data  and  Spectral  Analysis 

Figure  1  shows  a  representative  example  of  the 
temperatures  over  one  of  the  23  nights  of  observations.  The 
obvious  gaps,  due  to  realignment  procedures,  were  present  on 
all  nights  and  they  made  direct  spectral  analysis  impossible. 
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On  average  there  were  4  such  gaps  per  night  constituting,  on 
average,  about  1/4  to  1/5  the  duration  of  each  night's 
observation  interval.  Another  problem  was  that  the  sampling 
intervals  varied  somewhat  over  the  night.  This  however  was 
solved  by  cubic  interpolation  to  create  evenly  spaced  samples. 
Five  nights  had  ;ui  average  sampling  rate  of  19.9  sec.  (which 
differed  from  night  to  night,  o  =  0.2)  and  the  average 
observation  time  was  about  5  ±  1  hr.s.;  the  remaining  18  nights 
had  an  average  sampling  interval  of  39.5  sec.  (+.5)  and 
duration  4+1.5  hrs. 

Two  independent  methods  were  used  to  obtain  the  PSD's, 
The  first  was  the  well  known  one  of  Blackman  and  Tukey 
( 1959),  (BT),  and  in  this  case  the  problem  of  gaps  was  solved 
due  to  the  fact  that  there  were  no  gaps  in  the  autocorrelations. 
Ninety  jsercent  of  possible  lags  were  used  (instead  of  the  usu;il 
10%).  The  autocorrelations  were  obtained  by  the  zero-padded 
FFT  method  of  Oppenlieim  and  Schafer  (1975)  and  we  used 
the  "biased"  version  (normalized  by  total  number  of  data 
points)  to  avoid  the  subsequent  occurrence  of  negative  PSD's. 
Prior  to  autocorrelation  the  data  were  demeaned  and 
prewhitened  by  x.=x.-ox._|.  The  autocorrelations  were 
tapered  by  a  triangular  window  (Bartlett)  and  the  (whitened) 
PSD  obta'mei)  via  zero-padded  FFT.  "Post-darkening"  was 
obtained  by  dividing  tire  white  spectrum  by  ( 1  +  a- 
-  2a  cos  (coAt))  where  At  is  the  sampling  interval  (c.f. 
Blackman  and  Tukey,  (1959)). 

Straight  lines  were  least  square  fitted  to  both  the  white  and 
post  darkened  spectra.  The  values  of  a  =  1.0,  0.95,  and  0.90 
gave  whitened  .slopes  of  +0,3,  +0.006  and  -0,3  respectively 
hence  a  =  0.95  was  chosen  for  the  subsequent  analysis.  The 
fitted  lines  extended  from  a  minimum  frequency  (with  period 
equal  to  1/3  the  observation  interval)  to  the  maximum 
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frequency  of  2.5  x  10“-^  HZ  which  is  approximately  HO'IE  of 
the  Brunt  Vaisalla  frequency  (i.e.  maximum  gravity  wave 
frequency).  This  seemed  to  avoid  noise,  acoustic  compression 
and  other  contaminants. 

It  should  he  mentioned  that  we  have  made  the  usual 
a.ssumption  (made  in  lidar  and  ratku-  ohvservations)  that  witliin 
such  a  frequency  range  the  oscillations  are  due  exclusively  to 
gravity  waves. 

The  other  PSD  tecluiique  useii  was  the  maximum  entropy 
method  (MEM)  of  Burg  using  the  "Burg  algorithm"  (see 
Marple,  1987,  p.  213).  Tliis  was  applied  to  each  "piece”  of 
data  between  the  gaps  which,  in  effect,  supplied  an  alternate 
solution  to  the  gap  prohlem.  We  used  7  linear  predictor 
coefficients  the  choice  being  based  on  the  relatively  .slow 
decrease  of  prediction  error,  beyond  that  point. 

As  will  he  shown,  the  two  methods  gave  very  similar 
results  and.  in  addition,  they  both  gave  appropriate  responses 
to  a  family  of  known  synthetic  signals.  Figure  2  shows  a 
superposition  of  PSD's  ohiauied  both  methods  for  the  data 
in  Figure  1.  In  all  23  cases  the  MEM  result  looked  like  a 
smoothed  version  of  the  BT  result.  Straight  lines  were  fit  over 
the  same  ranges  in  the  MEM  PSD's  as  in  the  BT  PSD’s. 
Results  from  the  two  methods  did  not  perfectly  coinckle  on  an 
individual  basis,  but  there  was  very  good  overall  statistical 
agreement  as  will  be  seen. 

Theory  for  (w),  the  PSD  of  Time  Variations  of 
Temperature 

The  following  assumes  that  temperature  variations  between 
frequencies  (O  =  cOj  (inertial)  and  N  (buoyancy  frequency)  are 


:irz) 


Figure  2.  PSD’s  of  Temperatures  Given  in  Figure  I.  Both 
Blackman  and  Tukey  and  Burg  techniques  (smoother  one) 
were  used.  This  shows  an  average  case. 


due  to  gravity  waves  (mean  winds  are  ignoretll.  See 
Appendix  A  for  the  justification.  It  was  shown  in  Makhlouf  et 
al.  (1990)  that  two  separate  physical  mechanisms  account  for 
wave  intluced  temperature  variations.  One  of  these  is 
essentially  acoustic  in  nature  and  the  other  is  due  to  altitude 
variations  of  the  fluid  ptucels.  In  the  latter  case  the  change  in 
temperature,  AT.  due  to  a  change  in  altituile,  Az.  is  given  by 


where  H  is  the  scide  height  (which  is  5.7  km  at  85  km,  U.S. 
Standiird  Atm.,  (1976))  rurd  y  is  1.4  (specific  heat  ratio).  F-or 
our  purposes  Az  =  lOOm  would  give  us  approximately  AT  = 
1‘^K.  It  can  be  shown  that  the  above  mentioned  acoustic 
component  can  be  ignored  when  w  <  N/2  at  85  km  altitude.  At 
higher  frequencies  the  assumption  needs  further  justification. 
We  found  to  <  .8N  to  be  satsfactory  as  was  mentioned  above. 

How  can  the  above  be  employed  to  derive  'Ey  (to)  from 
theory’.^  In  Dewan  (1990,  1991 )  a  similitude  argument  based 
on  a  wave  cascade  hypothesis  led  to  the  frequency  PSD  for 
vertical  particle  velocities  due  to  gravity  waves.  4^^.  (to),  given 
by 

0  0  0 

r-  ,  ’  :  2  2 ' 

a  e  to  -  to 

'*'v  i . ,  (3,2) 

to"  I  !  N" 

ui  units  of  (m2/s2)/(Hz).  where  e  is  the  energy  dissipation  rate 
and  a'  is  a  constant  of  order  unity.  Term  (1)  is  the  wave 
energy  density  and  is  in  some  seu.,c  analogous  to  a  similar 
relation  given  in  Tennekes  and  Lumley  (1972)  in  the  context  of 
turbulence.  Temi  (2)  is  the  fraction  of  wave  energy  in  the 
form  of  vertical  velocity  variance  at  to  and  is  derivable  from  the 
polarization  relations  (cf.  Garrett  and  Munk  (1975)  and 
Scheftler  and  Liu  ( 1985)).  The  third  tenn  converts  radians/sec 
to  Hz.  It  can  be  shown  that  (using  Az  =  1 00m  implies 

rK  =  AT); 


4*  (to) 

T.  (0))=  v  =4>  ( (0)71 00m 


(compare  Garrett  and  Munk  (1975)),  thus  for 
to  »  tO;  we  have 


a'27tf  7  "k 

4't.  to)  =  - 

T  ^,2:1  lOOm 

N  to 


Appropriate  values  for  the  pinameters  are  N  =  0,02  rad/sec 
(Dewan  ;tnd  Go..d  (1986)),  and  a'  =  2.2  (Dewan,  1990). 

Eq.  (3.4)  thus  relates  E  to  measured  values  of  4'^  (to)  ;uid  vice 
versa.  For  this  rea.son,  if  the  above  is  valid,  it  may  be  possible 
to  e.stimate  average  E  over  certain  altituile  regions  by  means  of 
mterferometric  measurements  of  airglow  temperature  flue- 
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tuations  (e.g.  OH  at  83  km,  Na  at  90  ktri,  ajid  01  at  O*'  km). 
Equally  promising  me  the  new  lidar  temperature  nteasurement 
techniques  described  by  Von  Zahn  (1987),  Gardner  (1991), 
and  She  (1991).  In  the  meantune  it  is  essential  to  test  the 
theory  by  experimentally  measuring  6  at  the  same  time  that 
Ty  (0))  is  obtained.  One  way  to  do  this  is  by  means  of 
Doppler  broadenitig  of  radar  returns  as  described  in  Hocking 
(1988).  Radar  measunnents  of  4*^,  (to)  would  be  very  useful 
in  addition  of  course. 


Comparison  Between  Tlieory  ami  Exireriment 


Figure  3  summarizes  all  the  fitted  lines  (23)  for  the  BT- 
method  derived  V-j-  (co)  spectra.  Tlie  average  slope  was  -1.7  + 
0.3.  The  MEM  spectra,  while  very  .smooth,  gave  similar  line 
fits  with  slopes  -1.7  +  0.2  on  average.  These  results  are 
consistent  with  the  theoretical  value  of  -2.0.  (On  the  other 
hand  if  a  larger  data  set  were  to  narrow  the  spread  around  - 
1.7,  the  proverbial  -3/3  slope  would  be  a  more  valid  choice 
(cf  the  appendix  of  Dewan,  1990)).  We  now  compare  V-j-  (co) 
amplitude  with  the  theoretical  connection  with  observed  £ 
values. 

Solving  Eq.  (3.4)  for  e. 


£  =  (CO)  (0 


2 


N‘  (100)'’ 

2n  a' 


(4.1) 


=  (0.29) ;  iir.|.  (co)  •  co'j 

Figure  4  shows  the  23  fitted  lines  of  M'j  (co)  •  C0“,  As  can 
be  seen  they  fall  mainly  in  the  amplitude  range  between  10"- 
and  10“l  which,  from  Eq.  (4.1)  gives  the  approximate  range 
(0.003  <  £  <  0.03  m-/s-^).  This  overlaps  some  observations  of 
£  (at  another  time  and  geographic  location)  by  Hocking  (1988) 
in  the  altitude  range  of  80-90  km.  These  results  gave  (0.01  < 
£  <  0.1  m^/,s3).  This  outcome  is  sufficiently  encouraging  to 


Figure  3.  Summary  of  all  23  PSD’s.  Fitted  lines  from  BT 
derived  PSD's  used  here. 


Figure  4.  Sununary  of  BT  Fitted  Lines  for  Vy  (to)  •  to^. 

suggest  that  it  would  be  most  desirable  to  conduct  a  direct  test 
by  means  of  simultaneous  radar  (£)  and  interferometer  (or 
lidtn)  (T)  measurements. 

A  word  of  caution  is  needed.  Under  certain  circumstances 
the  measured  (co)  is  not  flat  (c.f.  Eq.  3.2)  but  instead  can 
have  a  negative  slope  (Ecklund  et  al.  (1986)  and  Sato  (1990)). 
Why  Ty  (to)  occasionally  has  a  negative  slope  (-3/3)  has  not 
yet  been  fully  explained.  Since  Eq.  (3.2)  holds  only  in  a  fluid 
at  rest  (i.e.  intrinsic  coordinates)  David  Fritts  (in  private  con¬ 
versation)  has  suggested  that  Doppler  effects  may  lie  the  cause 
of  the  observed  negative  slopes.  Needless  to  say  the  results  of 
the  present  paper  would  not  be  valid  under  those  conditions 
where  Eq.  (3.2)  is  not  valid.  See  also  Van  Zandt  et  al  (1991 ). 

Conclusions 

OH  rotational  temperatures  were  measured  (for  23  nights) 
and  the  'Ey  (co)  PSD's  obtained.  The  results  appear  to  agree 
with  the  theory  of  internal  gravity  wave  PSD’s  presented  in 
Dewan  (199d,  1991)  and  deduced  values  of  £  overlap  the 
range  experimentally  measuretl.  Direct  radar  and 
interferometer  tests  were  suggested  and  if  they  result  in 
confirmation  of  the  theory,  one  will  have  further  evidence  that 
the  wave  cascade  idea  first  suggested  in  Dewan  (1979)  is 
valid.  We  would  also  have  a  new  method  to  obtain  £  via 
interferometer  or  lidar  temperature  measurements  (under 
appropriate  conditions). 
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Appendix  A 


During  the  1988  Colorado  campaign  Prof,  S,  Avery  of  the 
LIniversity  of  Colorailo  made  meteor  wind  measurements  from 
the  Platteville  radar  site  which  was  centrally  located  with 
respect  to  our  instmments.  The  device,  known  as  MEDAC 
(Wang  et  al.  1987),  was  used.  In  this  way  we  obtained  a 
mean  wind  altitude  profile  for  days  127  to  142,  and  another 
average  for  199  to  205.  These  intervals  covered  16  of  the  23 
nights  of  observation.  At  no  time  did  these  me:m  winds,  u. 
exceed  15  m/s  in  the  important  OH  emitting  range  of  83-87 
km.  We  now  assume  this  to  lie  at  least  approximately  valid  for 
the  entire  campaign. 

According  to  Van  Zandt  el  al  (1991)  the  parameter 
u.  2  u/c,  (where  c,  is  the  characteristic  horizontal  phase  trace 
velocity)  detennines  the  doppler  effect  on  temporal  spectra. 

llsing  c.  =  ^  and  k^  =  from  Smith  et  al  (1987) 

‘  20  X  10  m 

one  obtains  u^  <  0.25.  Figure  8  of  Van  Zandt  et  al  (1991) 

shows  that,  in  this  case,  there  would  lie  negligible  effects  on 
our  spectral  measurements, 
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